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Abstract We present GalevNB (Galev for A-body simulations), an utility that converts 
fundamental stellar properties of A-body simulations into observational properties us¬ 
ing the GALEV (GAlaxy Evolutionary synthesis models) package, and thus allowing di¬ 
rect comparisons between observations and A-body simulations. It works by convert¬ 
ing fundamental stellar properties, such as stellar mass, temperature, luminosity and 
metallicity into observational magnitudes for a variety of filters of mainstream instru¬ 
ments/telescopes, such as HST, ESO, SDSS, 2MASS, etc.), and into spectra that spans 
from far-UV (90 A) to near-IR (160 pm). As an application, we use GalevNB to inves¬ 
tigate the secular evolution of spectral energy distribution (SED) and color-magnitude 
diagram (CMD) of a simulated star cluster over a few hundred million years. With the 
results given by GalevNB we discover an UV-excess in the SED of the cluster over the 
whole simulation time. We also identify four candidates that contribute to the EUV peak, 
core helium burning stars, thermal pulsing asymptotic giant branch (TPAGB) stars, white 
dwarfs and naked helium stars. 

Key words: stars: kinematics and dynamics—stars: C-M diagrams—stars: AGB and 
post-AGB —stars:white dwarf 

1 INTRODUCTION 

Models of dense stellar systems, such as globular or young dense stellar clusters, or nuclear star clusters 
with or without massive central black hole, require direct A-body simulations, which at least resolve 
and follow all stellar orbits in a star cluster as precisely as possible given the usable hardware and com¬ 
puting time. But due to the strong dependence on initial conditions of the gravitational A-body prob¬ 
lem, which originates from close encounters, the system is physically and numerically unstable (Miller 
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1966). Therefore results of direct high-accuracy A^-body simulations should always be carefully ana¬ 
lyzed, as they represent individual realizations of physically possible solutions. In another word, small 
variations of initial conditions may lead to different results with regard to individual objects or events. 
Nevertheless, the coarse-grained evolution (densities, distribution functions) is quite well following ex¬ 
pectations from statistical mechanics (see e.g. Giersz & Spurzem 1994, Spurzem 8c Aarseth 1996). And 
the simulations are an excellent tool to examine and study physical processes in star clusters. 

The sequence of high-accuracy direct A-body codes developed by Sverre Aarseth (Aarseth 1999; 
Aarseth 2003) is among the most widely used if not the most widely used code. Its most notable features 
are the Hermite 4* order integration scheme, the hierarchically blocked individual time steps (Makino & 
Aarseth 1992), the Kustaanheimo-Stiefel (KS) regularization for strong interactions (Kustaanheimo & 
Stiefel 1965) and its generalizations for few-body subsystems (chain regularization, Mikkola & Aarseth 
1993), and the Ahmad-Cohen neighbor scheme (Ahmad & Cohen 1973). NB0DY6++ is the current mas¬ 
sively parallel version of NB0DY6 (Spurzem 1999, Spurzem et al. 2008), and NB0DY6++GPU adds to it 
the use of accelerating many-core hardware (GPU) on every node (Wang et al. 2015). In the cited paper 
it becomes clear that direct A-body simulations with a million particles are feasible now. The challenges 
of the resulting massive simulation data and their possible approaches of management are addressed in 
|Cai et al. (20l4)l 

However, the output parameters of N-BODY simulations are mostly theoretical values. To make 
a direct comparison between our simulation data and observations, we combine GALEV (GALaxy 
Evolutionary synthesis models; Kotulla et al. 2009), a flexible algorithm to combine astrophysical col¬ 
ors in many Alters and spectra of stars (Lejeune, Cuisinier & Buser 1997, 1998) or sets of stars, with 
NB0DY6++GPU simulations. In this paper, we present this new code: GalevNB (Galev for A-body simula¬ 
tion) and its application in NB0DY6++GPU simulation data. Adapting subroutines from GALEV, GalevNB 
can produce spectra spans the range from far-UV (FUV) at 90 A to far IR at 160 pm, with a spectral 
resolution of 20 A in the UV-optical and 50-100 A in the near IR range. Given a list of requested Alters 
in HST, ESO, SDSS, 2MASS etc., GalevNB convolves the spectra with the Alter response functions and 
applies the chosen zero-points (Vegamag, ABmag, and STmag) to yield absolute magnitudes. GalevNB 
bridges theoretical parameters and their observed values, thus allows us to understand the color and 
spectra evolution of star clusters, and to determine the initial conditions and parameters of star cluster 
simulations with a direct comparison to observations. 

Though most Galactic and many Local Group clusters are resolved, observing individual stars in 
extragalactic star clusters is extremely challenging due to severe crowding. Therefore, integrated pho¬ 
tometry and spectroscopy have been used to identify extragalactic star clusters and obtain parameters 
(Bica et al. 1996a, Sarajedini et al. 2007, Peng et al. 2008, 2009, Johnson et al. 2012), and even to in¬ 
vestigate multiple stellar populations (Peacock et al. 2013). Using the imaging obtained with the High 
Resolution Channel of the Advanced Camera for Surveys on board HST, Larsen et al. (2011) man¬ 
age to construct color magnitude diagrams (CMDs) for several extragalactic star clusters within 5 Mpc. 
However, for the very distant star clusters, such as the ones in Antennae galaxies (~ 22 Mpc, Whitmore 
et al. 2007, Bastian et al. 2009), or in Virgo cluster (~ 16.5 Mpc, Peng et al. 2008, Peng et al. 2009), 
studies of individual stars are still not possible. Integrated colors can be used to derive the age of star 
clusters generally (Elson & Fall 1985). However, metal variation, stochastic effect (Girardi et al. 1995), 
and dynamical evolution (Fleck et al. 2006) will change integrated colors. Therefore, integrated pho¬ 
tometry may mot be enough to decouple the stellar and dynamical effects of distant star clusters. With 
GalevNB producing observational magnitudes and spectra for A-BODY simulations, it allows us to in¬ 
vestigate stellar evolution and dynamics (via colors and spectra) in distant star clusters at the same time, 
and even make prediction for observations. 

In this paper, we carry out star cluster simulations with NB0DY6++GPU codes, and use GalevNB 
to produce observational data for the simulated cluster. The computational methods of NB0DY6 +-hGPU 
simulations are summarized in Section An introduction to the structure and execution procedure 
of GalevNB is presented in Section |3] The observational features of the simulated cluster, CMDs and 
spectral energy distributions (SEDs) are outlined in Section |4] Finally, we present our discussion and 
summary in Section|3 
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2 COMPUTATIONAL METHOD AND INITIAL CONDITIONS 

The model clusters in this work are evolved using NB0DY6++GPU which is the MPI parallel version 
based on the state-of-the-art direct A-body integrator NB0DY6GPU (Aarseth 2003, Nitadori & Aarseth 
2012). Gravitational forces are computed using a fourth-order Hermite integration scheme with block 
time steps (and without softening). The code includes algorithms for stellar and binary evolution (Tout et 
al. 1997, Hurley et al.2000) and deals directly with perturbations to binary orbits, collisions and mergers, 
formation of three- and four-body subsystems, exchange interactions and tidal capture. The treatment 
of close encounters constitutes a large part of the code. 

We set up a simple dynamical model of a massive stellar cluster, as an example for GalevNB to work 
on. The stellar system is initially gas-free, and is assumed in virial equilibrium as an approximation 
after losing gas (i.e. the ratio of kinetic to potential energy is Qvir = 0.5). The single-aged population 
is run with Nq = 10000 single stars. Million-particle simulations of star clusters will be shown in 
our coming paper. The IMF is set up according to Kroupa (2001) with fixed lower and upper mass 
limits, m\ =0.1 Mq and mu = 20 Mq, respectively, resulting in an initial cluster mass of Mq - 4.7 x 
10^ Mq. To simplify the computation, stars are distributed according to Plummer model with a scale 
radius of 0.59 pc, which does not differ significantly from King model with Wo = 6, except in the 
outermost regions with low density. The models have an initial half-mass radius /?hm = 0.76 pc, and 
is not primordially mass segregated. We set the metallicity to sub-solar metallicity Zq = 0.001, which 
is similar to the halo population of Galactic globular clusters. Note that our simulated star cluster is 
still less massive than a typical Galactic globular cluster. But here we are interested in special spectral 
features, such as the production of UV bright stars, and aim to relate the current simulations to future 
million particle simulations. 

Since our primary interest is the short-term stellar evolution, we therefore carry out a simulation 
of 3000 A-body time units (corresponding to 655.6 Myr) with dynamical and internal stellar evolution. 
Snapshots of the simulation are selected to display in this paper for the presence of UV-bright stars (see 
Section 4). 

3 GALEVNB: GALEV FOR A-BODY SIMULATIONS 

3.1 GalevNB Structure 

In this section, we introduce the GalevNB structure and its parameters. The main program of 
GalevNB is GalevNB. £90, which parses single snapshot files (stellar evolution only) generated by 
NB0DY6++GPU (NB0DY6++)/ NB0DY6GPU (NB0DY6). It invokes seven subroutines (startomaginit, 
specint.initialize, reset.weights, startomag, add.star, spec2mag, spec.output)of 
GALEV package to convert effective temperature, stellar luminosity, metallicity, and mass into observa¬ 
tional magnitudes and spectra. The functions of these routines are presented in Table 1. The GalevNB 
package contains four folders: 1) spectral.templates, in which all the spectral template files from 
the BaSeL library of model atmospheres (Lejeune, Cuisinier & Buser 1997, 1998) are stored; 2) 
standard_filters, containing a large set of filter response functions (FUV, NUV, U, B, V, R, I, J, H, 
K) that are used as standard reference filters; 3) filter_response_curves, including filter response 
functions from magnitude systems of HST, ESO instruments, 2MASS, SDSS, Johnson, and Cousins 
in separate subfolders. We also provide a choice of user-specify filter response functions. Information 
about the entire set of available filters is presented in the file filterlist.dat. Please aware that 
filterlist. dat, in which the user specify their own choice of magnitude system by uncommenting 
the row of chosen filter, MUST be presented in the same directory as the NB0DY6++GPU (NB0DY6++)/ 
NB0DY6GPU (NB0DY6) snapshot files. The content of the file, filterlist. dat, is presented in Table 2. 

3.2 Installation and Usage of GalevNB 

To compile GalevNB, the user should have C-i-i- and fortran installed. The input file of GalevNB should be 
a sinlge snapshot output from NB0DY6++GPU (NB0DY6++)/ NB0DY6GPU (NB0DY6) simulations. In case of 
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Table 1 Functions of subroutines converting theoretical parameters into observational mag¬ 
nitudes and spectra 


Subroutine 

Function 

specint.initialize 

initialize the stellar spectra 

reset.weights 

reset the weight of stellar spectra 

add.star 

integrate the flux of all stars in the cluster 

spec_output 

output spectra 

startomaginit 

initialize the stellar magnitude 

spec2mag 

convolve the stellar spectra with the filter response function 

startomag 

compute magnitudes for stars 


Table 2 Column contents for the filter information file; filterlist. dat 

Column Content ID of zero point 

1 Filter name 

2 Corresponding path of the filter response function 

3 ID of selected zero point (default value is 1) 

4 Standard zero point in the Vega magnitude system 1 

5 Standard zero point in the AB magnitude system 2 

6 Standard zero point in the ST magnitude system 3 

7 Optional user-defined zero point 4 


a file containing all snapshots (called sev. 8 3 in NB0DY6++GPU/ NBODY6-t-+ and f ort. 8 3 in NB0DY6GPU 
/ NB0DY6), we provide the user with a shell script generate_snapshots. sh in the folder, scripts, 
for retrieving single snapshot data out of sev. 83 and fort .83. The user can select his/her preferred 
filters (maximum 20) by uncommenting the row of the corresponding filter in filterlist.dat, and 
choose his/her desired magnitude system (Table 2). Magnitudes of individual stars and the whole cluster, 
and spectra of the cluster or chosen stellar types are produced, respectively. The code of GalevNB is 
published onlineQ. Users can download it through internet. We also provide examples of output files 
from NB0DY6++GPU and GalevNB on the wet0. 

4 OBSERVATIONAL FEATURES OF SIMULATED STAR CLUSTERS 

With NB0DY6++GPU simulations done in Section 2, and the computation of GalevNB, we are able to 
investigate early evolution of CMD and SED of a stellar population in a real dynamical environment. The 
output of NB0DY6++GPU simulations contain integers to represent stellar type information (Hurley et al. 
2000) of each individual stars (Table 3). This allows one to identify exotic objects in the star cluster based 
on the evolution of CMD and SED. It is true that cluster dynamics, such as close stellar encounters or the 
formation of binaries, can significantly alter the stellar populations in dense stellar systems (Djorgovski 
& Piotto 1993). Nevertheless, to simplify the NB0DY6++GPU simulations that GalevNB works on, we do 
not consider primordial binaries in the current simulations. Therefore, stellar type 19-22 (binaries) do 
not present in our simulations. Due to the faintness of stellar remnant, stellar types 13-15 (though appear 
in our examples) are not included in the studies of CMD and SED. 

4.1 SED & CMD evolution 

After integrating fluxes of individual stars, GalevNB output total SED of the star cluster. In the main 
program, GalevNB. £90, user can choose to output SED of certain stellar types. Therefore, specific stars’ 
contribution to the total SED of the cluster can be quantified. In Eigures 1, 3, 5, 7, we display SEDs of 

* http://siIkroad.bao.ac.cn/repos/galevnb 
^ http://silkroad.bao.ac.cn/" xiaoying/ 
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Table 3 Stellar type defined in the NB0DY6++GPU codes 


Integer 

Stellar type 

0 

main sequence (MS) stars with mass< 0.7 MO 

1 

MS with mass> 0.7 MO 

2 

Hertzsprung Gap 

3 

Red Giant 

4 

Core Helium Burning 

5 

First (early) Asymptotic Giant Branch 

6 

Second Asymptotic Giant Branch (~Thermal Pulsing AGB) 

7 

Naked Helium Star MS 

8 

Naked Helium Star Hertzsprung Gap 

9 

Naked Helium Star Giant Branch 

10 

Helium white dwarf 

11 

Carbon/Oxygen White Dwarf 

12 

Oxygen/Neon White Dwarf 

13 

Neutron Star 

14 

Black Hole 

15 

Massless supernova remnant 

19 

Circularizing binary (c.m. value) 

20 

Circularized binary 

21 

First Roche stage (inactive) 

22 

Second Roche stage 


the cluster (black line), of asymptotic giant branch (AGB) stars (red line), of naked Helium stars (green 
line), of white dwarfs (blue), and of core helium burning stars (cyan line) at different ages. Movies of 
the evolution of SED and CMD are available in the weblink provided in Section 3.2^. 

At very early age, massive stars first exhaust their hydrogen in the core and ignite core helium 
burning. These core helium burning stars are UV-bright (Figure 1) and very blue (Figure 2), which are 
also termed as blue loops (Stothers & Chin 1991, El Fid 1995). The UV peak in the SED of the cluster 
(< 1000 A) is mainly due to the presence of these stars (Figure 1). 

As the cluster evolves, the intensity of the UV-excess drops (Figure 3, 5, 7). Besides core helium 
burning stars, we find another three stellar types producing UV radiation, second AGB stars, white 
dwarfs, and naked helium stars. 

1) Early AGB stars (stellar type = 5) are very red, and therefore mainly radiate in the red filters 
(Girardi 2000, Marigo et al. 2008, Salaris et al. 2014). After the third dredge-up, helium shells flash 
repeats many times. At this stage, early AGB becomes thermal pulsing AGB (TPAGB; Marigo et al. 
2008, Hurley et al. 2000). In the NB0DY6++GPU code, the TPAGB are called ’’second AGB” (stellar type 
= 6). The radii of the TPAGB grow so large that mass-loss is significant. At this time, the TPAGB reach 
very high temperature (> 10“^ K) and begin to be luminous in the UV filters (Figure 3). Mass-loss will 
eventually remove all the envelopes of the TPAGB, which may be seen as planetary nebulae (see the 
CMD in Figure 4). 

2) White dwarfs are born after the TPAGB turning into planetary nebulae. Though white dwarfs 
cool down eventually (Mestel 1952, Fiebert 1980, Hansen 8c Liebert 2003), young white dwarfs are 
very hot and blue (Mestel 1952, Rauch et al. 2014, Torres et al. 2014), whose radiation is largely at the 
UV (< 2000 A, Figures 5 & 7). Even though their luminosity at the UV is not high, their contribution to 
the UV-excess of the cluster is long-term and might not be neglected, since stars evolve to white dwarfs 
continually. 

3) One short term radiator at the UV is the naked helium star, with extremely blue color over the 
main-sequence turn-off (Figures 7 & 8). They are massive stars losing hydrogen envelopes, the so- 
called Wolf-Rayet stars in observations (Crowther 2007, Shara et al. 2013). Since naked helium stars 
only appear at young age for a short time-scale, they cannot influence the UV-excess of the cluster in 
the long run. 
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5 DISCUSSION 

We present an integrated software solution, GalevNB, a translator that bridges NBODY (6++ / 6) sim¬ 
ulations (with or without GPU) and observations. We run a short-term dynamical model of a star cluster 
based on NB0DY6++GPU code, producing data for GalevNB. GalevNB computes observational magni¬ 
tudes and spectra of stars by theoretical parameters, stellar mass, temperature, luminosity and metallic- 
ity, which are generated by NB0DY6++GPU simulations. In the SED evolution of the simulated cluster, 
we found a UV-excess in the wavelength shorter than 2000 A. 

Similar phenomenon is found in elliptical galaxies and early-type spiral bulges whose SEDs sur¬ 
prisingly increase to shorter wavelengths over the range 2000 to 1200 A, (Code 1969, Code et al. 1972, 
Eaber 1983, Burstein et al. 1988, Kurucz 1991), in contrary to the expectation of an old stellar systems 
that are usually assumed to be dark in the EUV. This rise in the spectrum at wavelengths shorter than 
2000 A is called ”UV upturn” or ”UV-excess”. The UV-excess resembles the Rayleigh-Jeans tail of a 
thermal source with effective temperature larger than 20000 K (Hills 1971). During the last decades, 
many efforts have been made to find out ’’culprits” of the UV-excess. In contrast to our simulated clus¬ 
ter whose age is young (up to 655.6Myr), in most of the early-type galaxies young massive stars are 
absent (O’connell 1986, Welch 1982, Buzzoni et al. 2012). Therefore, old, hot, low-mass stars become 
a more popular choice. Extreme blue horizontal branch stars, with core helium burning and very thin 
envelopes, might be promising for explaining this phenomenon. These stars are found in both metal-rich 
and metal-poor star clusters (Rich et al. 1997, Eerarro 1998, Li et al. 2013, Buzzoni et al. 2012), which 
makes their origins somehow ambiguous. Theorists proposed another candidate producing a significant 
amount of UV radiation, post-AGB stars (Buzzoni & Gonzalez-Lopezlira 2008), a quarter of which un¬ 
dergo TPAGB phase (Lawlor & MacDonald 2006). Some studies suggested that new-born white dwarfs 
were hot enough to emit moderate amount of UV photons (Hills 1971, Bica et al. 1996b). Recently, 
based on detail computations of binaries, Han et al. (2007) conclude that the brown dwarf binary should 
be a more promising candidate accounting for the UV-excess. 

Most of the candidate stars mentioned above originate in the cluster environment. With the applica¬ 
tion of GalevNB to NB0DY6++GPU simulation data, we are able to observe the SED evolution of different 
stellar types. Through this way, we find out four kinds of hot, blue stars that are dominant in the UV- 
excess: core helium burning stars, second AGB (TPAGB), white dwarfs and naked helium stars. Among 
them, second AGB is a favorable candidate from theoretical point of view (O’Connell 1999). On the 
contrary, white dwarf’s candidate position is controversial (Magris & Bruzual 1993, Landsman et al. 
1998) because of low luminosity. The life time of naked helium stars is very short. Though they are very 
bright at the UV, their short-term emission makes them insignificant candidates. However, how the UV 
radiation of the candidate stars evolve in the long-term is beyond the scope of this paper, which cannot 
be achieved with a simple stellar model. A detailed NB0DY6++GPU simulation investigating the evolu¬ 
tion of UV-excess in star clusters, with more particles, realistic initial condition and primordial binaries, 
will be presented in our coming paper (Pang et al. in prep.). The NB0DY6++GPU code in its latest GPU 
accelerated version is publicly available in Wang et al. 2015. We provide our scripts and subroutines 
of GalevNB online'. Please notice that these are compatible with both NB0DY6++ & 6 (with or without 
GPU) version. Users of the other version of NBODY may have to adjust some data formats. 

Eor more than a decade ’’simulating observations of simulations” has been a topic in the community 
of MODEST (MODelling Dense Star Clusters, see www.manybody. org/modest/ and Kouwenhoven 
et al. 2004). Hurley et al. (2005) have published a pioneering study, which presented a full Hertzsprung- 
Russell diagram (luminosity vs. effective temperature) of all stars in the cluster M67 obtained from a 
direct A-body simulation with NB0DY6, to be directly comparable to observations. These capabilities 
come with the standard public versions of NB0DY6 and NB0DY6++ (with or without GPU) already. The 
Monte Carlo code MOCCA is now also providing similar capabilities (Giersz, Heggie & Hurley 2008); 
recently it has been extended to allow simulated observations of star clusters with specific telescopes 
using their colour windows and standard software used by observers (Askar et al. 2015). This work has 
been driven by future key observational projects (for example by Hubble Space Telescope) for globular 
clusters (e.g. Milone et al. 2014); however, to fully uncover the 6D dynamical structure and the 3D 
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gravitational potential, as well as the population history of star clusters, pure photometry may not be 
sufficient. Future deep high resolutions spectroscopic and spectrophotometric observations of globular 
and other star clusters are required in order to uncover one more dimension in velocity (radial velocity) 
and obtain more information about the stellar population through full spectra. Kacharov et al. (2014) and 
den Brok et al. (2014) provide examples on how kinematic data provide key insights into the dynamics 
of globular clusters (rotation and the quest for intermediate mass black holes). See also for an overview 
van de Ven (2010). Here the combination of the GALEV codes, previously mainly used for synthesis of 
galactic stellar populations (Kotulla et al. 2009) with our most recent direct A7-body code for star clusters 
is a pioneering step to provide the corresponding full spectral information at any time and any region in 
the cluster, star-by-star or in integrated fields from the models. 
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X. Y. Pang & C. Olczak & D. F. Guo & R. Spurzem 



Fig. 1 SED of the simulated cluster (black line) and core helium buring stars (cyan line) at 
the age of 14.4 Myr. SED of other stellar types will be shown in the following figures. SED of 
core helium burning stars (stellar type = 4) is indicated as cyan line, AGB stars (stellar type 
= 5 & 6) red line, naked helium stars (stellar type = 7, 8, 9) green line, white dwarfs (stellar 
type = 10, 11, 12) blue line. High resolution figures are available for online version. 


Plummer Cluster CMD 14.4 Myr, N = 9990 
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Fig. 2 J-H Color-Magnitude diagram of the simulated cluster at the age of 14.4 Myr. N is the 
total number of stars, and n is the number of stars of each stellar type indicated in the separate 
box to the right. Stellar types (see Table 3) are marked with different color and by different 
symbols. Main sequence stars (Light blue filled round dot; stellar type = 0; sky blue one-pixel 
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Fig. 3 SED of the simulated cluster (black line), AGB stars (red line), core helium burning 
stars (cyan line) and white dwarfs (blue line) at the age of 216.6Myr. Color coding is the 
same as Figure 1. High resolution figures are available for online version. 


Plummer Cluster CMD 216.6 Myr, N = 9527 



Fig.4 J-H Color-Magnitude diagram of the simulated cluster at the age of 216.6Myr. N is 
the total number of stars, and n is the number of stars of each stellar type indicated in the 
separate box to the right. Stellar types (see Table 3) are marked with different color and by 
different symbols. Main sequence stars (Light blue filled round dot: stellar type = 0; sky blue 
one-pixel point: stellar type =1), core helium burning stars (light green left-pointing-triangle: 
stellar type = 4), AGB stars (green downward-tripple-point: stellar type =5; green upward- 
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Fig. 5 SED of the simulated cluster (black line), AGB stars (red line), core helium burning 
stars (cyan line) and white dwarfs (blue line) at the age of 466.8 Myr. Color coding is the 
same as Figure 1. High resolution figures are available for online version. 


Plummer Cluster CMD 466.8 Myr, N = 8651 



Fig. 6 J-H Color-Magnitude diagram of the simulated cluster at the age of 466.8 Myr. N is the 
total number of stars, and n is the number of stars of each stellar type indicated in the separate 
box to the right. Stellar types (see Table 3) are marked with different color and by different 
symbols. Main sequence stars (Light blue filled round dot: stellar type = 0; sky blue one-pixel 
point: stellar type =1), Hertzsprung gap stars (dark blue filled downward-triangle: stellar type 
=2), core helium burning stars (light green left-pointing-triangle: stellar type = 4), AGB stars 
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Fig. 7 SED of the simulated cluster (black line) naked helium stars (green line), core helium 
burning stars (cyan line) and white dwarfs (blue line) at the age of 504.4 Myr. Color coding is 
the same as Figure 1. High resolution figures are available for online version. 


Plummer Cluster CMD 504.4 Myr, N = 8493 



Fig. 8 J-H Color-Magnitude diagram of the simulated cluster at the age of 504.4 Myr. N is 
the total number of stars, and n is the number of stars of each stellar type indicated in the 
separate box to the right. Stellar types (see Table 3) are marked with different color and by 
different symbols. Main sequence stars (Light blue filled round dot: stellar type = 0; sky blue 
one-pixel point: stellar type =1), naked helium burning stars (pink-red square: stellar type = 
8), core helium burning stars (light green left-pointing-triangle: stellar type = 4), and white 
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